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We apply classical algorithms for approxi¬ 
mately solving constraint satisfaction problems 
to find bounds on extremal eigenvalues of lo¬ 
cal Hamiltonians. We consider spin Hamilto¬ 
nians for which we have an upper bound on 
the number of terms in which each spin par¬ 
ticipates, and find extensive bounds for the op¬ 
erator norm and ground-state energy of such 
Hamiltonians under this constraint. In each 
case the bound is achieved by a product state 
which can be found efficiently using a classical 
algorithm. 

1 Introduction 

The eigenvalue statistics of a local Hamiltonian are 
related to its structure. One example is the level 
spacings of chaotic vs integrable systems, which can 
be seen as the small-scale structure of the spectrum. 
What about the large-scale features, such as the ex¬ 
tremal eigenvalues? Do these scale differently for 
non-interacting or interacting systems? It is gener¬ 
ally understood that interacting systems can be frus¬ 
trated, meaning that all the local terms cannot simul¬ 
taneously be in their ground state. This situation is 
generically true for local terms with entangled ground 
states. But how much of an effect can frustration have 
on the ground-state energy of a system? 

Here we study the extremal eigenvalues of quan¬ 
tum Hamiltonians which are only weakly interacting, 
in the sense that they can be written as sums of terms 
where each term depends only on a few qubits, and 
each qubit is included in only a few terms. With this 
mild form of locality imposed, how far apart must the 
largest and smallest eigenvalues be? If the Hamilto¬ 
nian were non-interacting, the separation should scale 
with the size of the system. For a more general Hamil¬ 
tonian, the extremal eigenvectors may be highly en¬ 
tangled and interacting terms may contribute oppo¬ 
site signs. Nevertheless, in this paper we show lower 
bounds on the norms of local Hamiltonians under very 
general conditions. An additional argument shows 
specifically that the ground-state energy is low (or 
if desired, that the top eigenvalue is high). 
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Theorem 1. Let H be a traceless k-local Hamiltonian 
on n qubits such that k = 0(1). Assume that H can 
be expressed as a weighted sum of m distinct Pauli 
terms such that each term has operator norm 0(1), 
and each qubit participates in at most £ terms. Then 
||iJ|| > and XminiH) < In each 

case the bound is achieved by a product state which 
can be found efficiently using a classical algorithm. 

In the above theorem, ||il|| is the operator norm of 
H and is the lowest eigenvalue (ground-state 

energy) of H. (Of course a similar statement could 
also be made about Amax- We focus on Amin because of 
its relevance to physical systems and to constraint sat¬ 
isfaction problems.) The notation 0{f{x)),fl{f{x)) 
refers to functions that are < cf{x) or > cf{x), re¬ 
spectively, for some absolute constant c. We write 
0(/(x)) to mean both a function that is both 0{f{x)) 
and Ll{f{x)). 

If H were a non-interacting Hamiltonian {k = £ = 
1) the largest and smallest eigenvalues would both be 
0(m) = 0(n). Thus Theorem 1 can be viewed as 
saying that interaction can reduce the norm of H by 
at most a O(v^) factor and can reduce the smallest 
eigenvalue by at most a 0{£) factor. Observe that 
the bound on Xmin{H) is —Lt{n) for all lattice Hamil¬ 
tonians. This proves that for any such system the 
ground-state energy is smaller than the average en¬ 
ergy by an extensive amount. By constrast, using our 
information about tr alone would only show that 
Amin(?^) ^ —Ll{^/rn), which is in general a vanishing 
fraction of system size. 

The restriction to terms of weight 0(1) in Theo¬ 
rem 1 is not essential and is only included to simplify 
the bounds. Further, the hidden constants are not 
overly large for small k; more precise statements of 
our results are given below. For example, for 2-local 
qubit Hamiltonians, the precise bound on Amin we ob¬ 
tain is Amin(-f?) < ~\\H\\i/{‘24:£), where ||i?||i is the 
sum of the absolute values of the coefficients in the 
Pauli expansion of H. As a simple instance where this 
bound can be applied, consider the antiferromagnetic 
Heisenberg model H = XiXj + YiYj + ZiZj on 

a regular lattice with n vertices. Then, for any such 
lattice, we obtain Xmin{H)/n < —1/48. 

Theorem 1 can be applied to qudit Hamiltonians 
with local dimension d > 2 by embedding each sub¬ 
system in |"log 2 d~\ qubits at the expense of increasing 
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the locality from k to k [log 2 d\. 

Proof outline. Both results that make up 
Theorem 1 are based on the use of a correspon¬ 
dence between local quantum Hamiltonians and low- 
degree polynomials, which allows us to apply classi¬ 
cal approximation algorithms for constraint satisfac¬ 
tion problems. This correspondence uses a qubit 2- 
design [5, 13] to convert arbitrary qubit Hamiltonians 
to polynomials over boolean variables. 

The operator norm bound in Theorem 1 (stated 
more precisely as Lemma 3 below) is based on re¬ 
cent work of Barak et al. [2] which gives an efficient 
randomised algorithm for satisfying a relatively large 
fraction of a set of linear equations over F 2 . The 
bound on Amin (stated more precisely as Lemma 5 be¬ 
low) is based on analysing a natural greedy algorithm 
which is similar to a classical algorithm of Hastad [14]. 
Our results can be seen as generalising these two clas¬ 
sical algorithms to the quantum regime. 

Other related work. Bansal, Bravyi and Ter- 
hal [1] have previously shown that, for 2-local qubit 
Hamiltonians iJ on a planar graph with Pauli inter¬ 
actions of weight 0(1), Amin(^f) < —0(m). Simi¬ 
larly to our result, their proof uses a mapping be¬ 
tween quantum and classical Hamiltonians and proves 
the existence of a product state achieving a —0 (to) 
bound. However, the two results are not comparable; 
ours holds for non-planar graphs and A:-local Hamil¬ 
tonians for fc > 2, while theirs encompasses two-local 
Hamiltonians on planar graphs with vertices of arbi¬ 
trarily high degree. The quantum-classical mapping 
used is also different. Finally, the constants in our 
results are somewhat better (for example, they ob¬ 
tain Xjnin{H)/n < —1/135 for the antiferromagnetic 
Heisenberg model on a 2D triangular lattice). 

This work was motivated by [2] (whose main re¬ 
sult is presented in Section 3). Ref. [2] in turn was 
inspired by [7, 8], which gives a quantum algorithm 
for finding low-energy states of classical Hamiltonians. 
The relative performance of these different algorithms 
(ours/ [2] vs. [7, 8]) is in general unknown, and it is 
also open to determine the extent to which [7] can 
be generalised to finding low-energy states of local 
Hamiltonians. 

One other related work is [3], which showed that 
when k = 2 and the degree of the interaction graph 
is large, then product states can provide a good ap¬ 
proximation for any state, with respect to the metric 
given by averaging the trace distance over the pairs of 
systems acted on by the Hamiltonian. In particular 
this means they can approximate the minimum and 
maximum eigenvalues. Both our result and [3] yield 
similar error bounds (ours are somewhat tighter), but 
in this sense apply to incomparable settings: [3] show 
that product states nearly match the energy of some 
other state (e.g. the true ground state) with possibly 
unknown energy while our paper puts explicit bounds 
on the maximum and/or minimum energy. 


Another way to think about our work is as show¬ 
ing that interacting spins must nevertheless behave in 
some ways like noninteracting spins. In this picture, 
some vaguely related work is [10, 4], which show that 
under some conditions lattice systems have a density 
of states that is approximately Gaussian. These re¬ 
sults are incomparable to ours, even aside from the 
different assumptions, because we put bounds on the 
extremal eigenvalues while they study the density of 
states and/or the thermal states at nonzero tempera¬ 
ture. Theorem 4 of Ref. [11] also bounded the density 
of states of fc-local Hamiltonians under general condi¬ 
tions, but in the opposite direction: i.e. putting upper 
bounds on how many eigenvalues could have large ab¬ 
solute value. 

Why product states? Ground states of local 
Hamiltonians may be highly entangled [9]. But our 
bounds on j|77j| and XminiH) are achieved only with 
product states. One reason for this in the case of 
j|i7j| is that we are using random states, and prod¬ 
uct states have much larger fluctuations than generic 
entangled states. Indeed the variance of {tplHlip) for 
a random unit vector j-^) is only 0(m/2”). It is an 
interesting open question to find a distribution over 
entangled states that improves the constant factors in 
Theorem 1 that we achieve with product states. 

Fourier analysis of boolean functions. We will 
need some basic facts from classical Fourier analysis of 
boolean functions [12]. Any function / : {±1}" —t K 
can be written as 

/(®) = ( 1 ) 

SC[n] 

where xs ■= W^^gXi and [n] := {l,...,n}. This 
is known as the Fourier expansion of /. Parseval’s 
equality implies that 

Var(/) := E,[/(x)2] - E,[f{x)f = ^ hS)\ (2) 

where the expectation is taken over the uniform dis¬ 
tribution on {±1}". In addition, /(0) = Ea;[/(x)]. 
The influence of the /’th coordinate on / is dehned as 

lnf,{f)=J2fiSr- 

S3j 

2 The quantum-classical correspon¬ 
dence 

Let 77 be a /c-local Hamiltonian which has Pauli ex¬ 
pansion 

77 = ^ Hs Si ^ S2 ^ Sn 

sG{/,X,F.Z}" 

for some weights 77^ that we can view as a Fourier 
expansion of 77 analogous to that in (1). Define the 
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norms \\H\\p := for p > 1. In order to 

apply classical bounds to extremal eigenvalues of H, 
we observe that the action of a fc-local Hamiltonian 
on product states corresponds to a low-degree poly¬ 
nomial. Define the following set of states [13, 16]: 


3 Operator norm bounds 

We will use the following result of Barak et al. [2], 
which is a constructive version of a probabilistic 
bound previously shown by Dinur et al. [6]: 


1^^++) = ^ (^^S + v^lO) + 6-/^3-V3|l)) , 
\’<P-+) = Z\iIj++), IV'-h-) = X\tp++), IV'—) = Y\^++). 

These four states form a qubit 2-design; equivalently, 
a symmetric informationally-complete quantum mea¬ 
surement (SIC-POVM) on one qubit [13]. This mea¬ 
surement was studied in detail in [16]. Geometrically, 
the states describe a tetrahedron within the Bloch 
sphere [5]. 

Then define the functions Xs ■ K, for 

se{I,X,Y,Z},hy 

Xs{x) = {tpx\s\lpx)- 

These functions are pleasingly simple: one can verify 
that 

Xiix) = 1 , xx(x) = ^,XYix) = xzix) = 

( 3 ) 

Split each x G {±1}^" into n consecutive blocks of 
length 2, written as a: = and define 

the function /// : {il}^" —)■ M by 

fH{.x) = (V’2:(i)|---(V'a;(")|-ff|V'x(i))---|V'rr(")) 

= XI H,X8Ax‘^^^)Xs2i.x^^'^)---Xs^{x^'"'^). 

As each x G {±1}^” corresponds to a state 
|V’x{i)) ■ ■ ■ |V'k("))i we immediately have the 
bounds > max„,g{±i} 2 „ ///(x) and 

Amin(^l) < min 3 ,g{±i} 2 r. fnix). We will now 

proceed to show bounds on max^jg^^ipTi fnix) and 
minj-g^iip.. fnix) by viewing fnix) as a polynomial. 

As H is fc-local and each function Xs (s ^ I) is a. 
monomial of degree at most 2, fn is a polynomial of 
degree at most 2k. Because the Fourier expansion of 
each function Xs contains only one term, each term 
in H corresponds to exactly one term in the Fourier 
expansion of fn- Indeed 

iA(s) = H«3-I^I/2, (4) 

where s e {I,X,Y, Z}^ and |s| = |{* : Si ^ J}|. This 
corresponds to identifying {I, X,Y, Z}'^ with subsets 
of [2n] in the natural way. Thus by eqns. (2), (3) and 
(4) we have 

Var(/^)= X (5) 

Inf,(/^)= X ^'3-1^1. 

S,SjT^I 


Theorem 2 (Barak et al. [2]). There is a univer¬ 
sal eonstant C and a randomised algorithm such that 
the following holds. Let f : {±1}” —> K &e a poly¬ 
nomial with degree at most k such that Var(/) = 1. 
Let t > 1 and suppose that Infi(/) < C~^t~^ for all 
i G [n]. Then with high probability the algorithm out¬ 
puts X G {±1}" such that \f{x)\ > t. The algorithm 
runs in time poly(TO, n, exp(fc)), where m is the num¬ 
ber of nonzero monomials in f. 

Recent independent work of Hastad [15] describes 
an alternative, randomised algorithm achieving a sim¬ 
ilar bound. 

Given the quantum-classical correspondence dis¬ 
cussed in the previous section, we can now apply The¬ 
orem 2 to fn to prove the following result, which is 
one half of Theorem 1. 


Lemma 3. There is a universal constant D and a 
randomised classical algorithm such that the following 
holds. Let H be a traceless k-local Hamiltonian given 
as a weighted sum of m Pauli terms such that, for 
all j, InfjV/jy) < /max- Then with high probability 
the algorithm outputs a product state {ip) such that 
|('(/'|i?|V’)l ^ Zl“^||i?|| 2 /-\/-fmax- The running time of 
the algorithm is poly{m,n,exp(k)). 

Proof. First observe that if we simply pick x G 
{j_l} 2 n uniformly at random and consider the cor¬ 
responding product state lipx), 

Ex[{i^x\H\i;x)^] = Varifn) > ^ 

by (5). In addition (see e.g. [12, Theorem 9.24]), as 
fn is a degree-2fc polynomial, 

Pi'[l/ff(a')| > \/Var(///)] > exp{-0{k)). 

X 


Therefore, simply picking exp(0(fc)) random prod¬ 
uct states of the form \ipx) achieves \{'tpx\H\'f’x)\ > 
||iJ|i2/3^/^ with high probability. Let E be a univer¬ 


sal constant to be chosen later. If /max ^ E 


— k I 


H\\l 


-fci 


HWl /as desired. 


then |(V'a:|Al|V’a:)| > (Vsi?) 
taking D = \/3E. So assume henceforth that /max < 
E-'^\\H\\l Let = fn/V^arifn) so Var(/),) = 1. 
Then 


Inf, (/)/) = 


Inf,(/g) 

Ya.i{fH) 


< 3 


k /max 


l^lli 


Set 


^ \\Hh 


. /rrAd-v ■ J-n^ ■ ( 
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where C is the constant in Theorem 2 and we choose 
E large enough for the first inequality to hold. Then 
the algorithm of Theorem 2 outputs x G such 

that \fH{x)\ > £'“'"/ 2 ||H|| 2 /V/max- Rcnormalising 
again by multiplying by y^Va.r{J h) > 

= \fH{x)\ > {VsE)->^\\H\\yvT^, 

which completes the proof. 

Note that the algorithm does not need to know 
whether /max is large or not, since it can simply try 
both strategies and see which one results in the larger 
value of |(V^|i/|^)|. □ 

The first part of Theorem 1 is now immediate from 
Lemma 3. Let H he & fc-local Hamiltonian with k = 
0(1) such that H is a sum of m distinct Pauli terms, 
each of weight 0(1), with each qubit participating in 
I terms. Then \\H\\\ = 0(m), Infj(///) = 0(£). 

4 Bounds on extremal eigenvalues 

We now describe an algorithm for bounding extremal 
eigenvalues which is weaker, but holds for both the 
largest and smallest eigenvalues. Once again, the al¬ 
gorithm is based on applying the quantum-classical 
correspondence in Section 2 to a classical algorithm. 
We first describe the classical algorithm, which is a 
simple greedy approach to find large values taken by 
a low-degree polynomial on the boolean cube. 

Given / : {±1}" —> K define W := W(f) = 
1/(5') I. Assume that, for all i € [n], 

\{T C [n] : /(T) ^ Q and i G T}\ < i. 

Consider the following algorithm, based on ideas 
of [14] but somewhat simpler: 

1. Find S such that |/(5)| is maximal. 

2. Substitute values for Xi, i G S, such that /(0) 
increases by at least |/(5)|. 

3. Repeat until / is constant; call this constant /end- 

It is not obvious that step (2) can be achieved, because 
there might exist T C 5 such that f{T) ^ 0. Define a 
function fs by fs{x) = J2 tcs f{T)xT- For each T C 

5 such that T 0 and each a G {±1}, Ea;,a:s=a[a:T] = 
0- ^°^x,xs=sgn{f{S))ifsi^)] = /(0) + l/('S')l, and there 
must exist some y achieving fs{y) > /(0) + |/(5)|. 
Searching over at most 2^ different values x is suffi¬ 
cient to find y. 

Lemma 4. When the above algorithm terminates, 

fend> m + W/{2ke). 

Proof. Let fj be the new function produced at the 
j’th stage of the algorithm, with /o = /. Let Mj 
be the value of |/j_i(5)| corresponding to the set S 
chosen at stage j. Then /j(0) = /j_i(0) -I- Mj and 


W^(/i) > R^(/}-i) ~ 2k£Mj. The latter inequality is 
shown as follows. For each i G S, there are at most 
£ subsets T such that fj_i{T) yf 0 and i G T. So 
there are at most kt subsets T such that T n 5 yf 
0. For each such T, the substitution of values Xi, 
i G S, implies that fj{T) is set to 0, and for some 
other subset T', fj{T') = /j_i(T') ± /y_i(T). These 
are the only coefficients modified by the substitution 
process. Thus W{fj) can only decrease by at most 
2|/j_i(T)| < 2Mj for each T such that 5 n T y/ 0. 
Rearranging we have Mj > (kF(/y_i) — W{fj))/{2k£) 
for j > 1 and hence 


fend = m+Y.M, 
j 


> m + 


= m + 


^Em-i)-m) 

3 

W 

Wi 


as claimed. 


□ 


The following lemma is now essentially immediate. 

Lemma 5. There is a universal constant E and a de¬ 
terministic classical algorithm such that the following 
holds. Let El be a traceless k-local Hamiltonian which 
can be written as a weighted sum of m distinct Pauli 
terms such that each qubit participates in at most I 
terms. Then the algorithm outputs a product state 
lip) such that {ip\H\ip) > E~^\\H\\i/£. The algorithm 
runs in time poly(m, n, exp(fc)). 

Proof. Apply the algorithm of Lemma 4 to the degree- 
2k polynomial fn ■ {±1}^" —> K defined as in Sec¬ 
tion 2. We have kF(/yy) > 3“^/^||//||i, /iy(0) = 0. 
Hence the algorithm finds |'0) such that > 

3-*'/2||.H|1i/(4M) > E-'^\\H\\i/£ for a large enough 
universal constant E. □ 


Applying the same procedure to —H is of course 
sufficient to also find \ip) such that {ilj\H\tp) < 
—E~'‘\\H\\i/£. Observing that ||i/||i = 0(m) for a 
Hamiltonian H which is a sum of m distinct Pauli 
terms of weight 0(1) completes the proof of the sec¬ 
ond part of Theorem 1. 


5 Optimality 

Both bounds in Theorem 1 are tight, even for clas¬ 
sical Hamiltonians, as demonstrated by the following 
examples. The first is based on an example in [2]. 
Consider the 2-local Hamiltonian on n qubits 

H= 'y ] Hij ZiZj, 

ihL&E 


where E is the edges of an arbitrary r-regular undi¬ 
rected graph on n vertices, and each weight Hij € 
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{±1} is picked uniformly at random. Then m = rnl2, 
£ = r. For each fixed x G {0,1}", 

{x\H\x)= Y. 

is a sum of rn/2 uniformly distributed elements of 
{±1}. Via a standard Chernoff bound argument, 

Pj[\{x\H\x) \ > t] < 2e"™. 

Fixing t = 0(nyV) and taking a union bound over all 
X G {0,1}", \\H\\ = 0(ny/r) = Oimlsft) with high 
probability. 

Second, consider the 2-local Hamiltonian on n 
qubits 

H = Yza,. 

i<j 

Then iJ is a sum of to = 0(n^) Pauli terms of weight 
1, where each qubit participates in £ = 0(n) terms. 
We have 

/ irri \ (n-2|a;|)^-n 

{x\H\x) = Yi-^) -•, 

i<j 

so Amin(^^) > —n/2 = —0(to/£). 

It is an open question whether both the bounds on 
||i?|| and Aniin(^^) can be saturated at the same time. 
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